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ABSTRACT 


The validity of the Glauber eikonal approximation has 
beecneassecced for clastic scattering ofmucleons musing 
Oputca Wapotentials torvlead, «calcium and helium ativarious 
energies from 100 MeV to 1 GeV. 

Since the above nuclei are even-even, the possible 
interactions are only central and spin-orbit, therefore an 
optical potential in the form of a complex central and 
complex spin-orbit modified Gaussian was used, the parameters 
of which were chosen so as to reproduce experimental data as 
closely as possible. There are three possible observables 
from such a system, and so three are studied here -- cross- 
section, polarisation and a polarisation transfer coefficient. 

PirsmtSomiicwiLrs tus tudaveot 1lsakind s10 COnsideraalio possible 
independant observables, most others stopping short at just 
considering the cross-section. 

Anvextira multiplicative factor has been found tor the 
spin-flip amplitude whose presence is required formally, but 
whose effect upon the approximation is not found to be 
SrrireL Cant. 

The results are presented graphically since each 
reader will have his (or her) own interpretation of when the 
approximation "fails". An interestig trend which emerged 


from the study is that at 0.5 GeV the approximation fairs 
aimogteas well as»at 2 GeV, and that the expected result of 


decreasing angular range with increasing atomic weight was 


found to be less severe than expected. 
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CHAPTER I 


INTRODUCTION 


In 1958 R.J. Glauber gave a lecture on an eikonal 
approximation for high energy collisions at the theoretical 
physics summer conference in Colorado (GL59). The approxima- 
tion, now named after him, started to take shape a few years 
before; the earliest relevant papers appearing in 1947, one 
by Moliére, (MO47), the other by Serber (SE47). The paper by 
Moliére seems to be the first to clearly relate the eikonal 
phase to the classical action integral of the particle's 
trajectory, and the paper also gives a discussion of the re- 
lated Born and W.K.B. approximations. 

In 1954 J.B. Malenka (MA54) gave the expression for 
the amplitudes for scattering by spin-3 particles, the non- 
spin-flip amplitude agreeing with Moliére's result. 

Watson (WA53) put forward a multiple scattering formal- 
ism in 1953 in which scattering and absorption of particles 
were treated using the one and two particle nuclear densities, 
obtaining a formal solution for the scattering ampltiude. By 
employing the impulse approximation, a solution is found which 
is of the form of a multiple scattering solution, and by con- 
Sidering only on-shell scattering for a large number of nucle- 
ons, he relates the formalism to the optical model. 

Glauber dgidea Jot. of work)in the field during’ the 1950"s 
in particular showing (GL55) the additive property of the ei- 
konal phases paving the way for the multiple scattering formal- 
ism, and by the end of the decade brought all the separate re- 
sults together to form what is now called "Glauber Theory." 

The term Glauber Theory in the literature has a double meaning, 
the eikonal approximation or the multiple scattering theory, 


in this study the concern is with the former. 
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The eikonal approximation may be thought of as the 
result obtained by converting the sum over phase-shifts into 
an integral over an impact parameter, as shown by Razavy et 
al (RA74), and is therefore sometimes referred to as an inm- 
pact parameter representation for the scattering amplitude. 
Other equivalent ways of obtaining the result are linearising 
the Hamiltonian, as is used in Chapter II of this thesis, or 
making the Schroedinger equation into a first order equation 


by means of the approximations 


+ «1 cnt Kae <4 


as is done by Glauber (GL59). In the above 9 is the scat- 
tering angle, V is the potential strength, E is the 
energy of the incident particle, k is the wavevector of 
theminerdenteparticle jw ali S a.leneth of ,order of «the “spatial 
extent of the potential, d is the typical minimum distance 
over which the potential varies significantly. To satisfy 
these three inequalities, the potential must be smooth in 
some sense. A clear exposition of the relationship between 
the last two methods of derivation may be found in a paper 
by Kamal (KA72). 

The linearised Hamiltonian approach is dealt with in 
quite some depth by Osborn (0S70), where off-shell effects 
are studied systematically and the linearised theory is shown 
to be exactly unitary (a consequence of the hermitian nature of 
the linearised Hamiltonian). 


Another way of looking at the eikonal approximation 
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is to expand the exponential in the expression for the 
scattering amplitude 


£(0) = ix [ag(an)ts - eX(b)) » ap 


(2) 
to obtain a series in powers of xX, which corresponds term 


for term with the Born Series. The first order terms in 
Cachescries turn Out vo be identical, and Byron.eu al. (BY73) 
have investigated relationships between higher order terms. 

There are several other high energy approximations, 
and it is interesting to compare the accuracy of each to that 
of the Glauber eikonal approximation. 

From the above it is apparent that the Glauber eikonal 
approximation is more accurate than the Born approximation, 
although Byron et al. (BY73) have pointed out that for a 
linear combination of Yukawa potentials, the second Born 
approximation is the most accurate at higher energies. 

Razavy et al. (RA76) have looked at an alternative way 
of linearising the Hamiltonian, obtaining an expression for 
the amplitudes which, upon numerical testing, is more accu- 
rate than Glauber's result, although the expression is com- 
plicated. 

Blankenbecler and Goldberger (BL61) have developed 
an impact parameter approximation.. The expression they have 
obtained for the amplitude, which agrees with the Glauber 
eikonal amplitude to order 56s is discussed and the formalism 
extended to spin-3 scattering by Kamal (KA72). 


T, Adachi et al. (AD65) have developed an impact 
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parameter formalism without approximation, but it yields 
a rather less manageable expression for the scattering 
amplitude. 

R; Sugar and R. Blankenbecler (SU69), Schiff (SC56) 
and Saxon and Schiff (SA57) have also developed eikonal for- 
malisms, although to date no one has found expressions for 
the amplitudes of comparable simplicity to those of the 
Glauber eikonal, which gives significantly better results. 

Y. Hahn (HA69)(HA70)(HA73) has carried out a survey 
of impact parameter formalisms, and comes to the conclusion 
that all fail beyond a certain value of momentum transfer. 
The approximation of Schiff (SC56) turned out to be the most 
accurate numerically, but the differences were slight and the 
survey was taken at only one energy and one geometry for the 
potential. 

There have been several attempts in the past to give 
corrections to the Glauber eikonal approximation, the most 
Bheren linn and successful approach seems to be that of Wal- 
lace (WA70)(WA71), who writes the amplitude in an expansion 
in powers of the difference between the exact and eikonal 
propagators, thereby obtaining a series with the Glauber 
eikonal approximation as the first term, which converges 
reasonably rapidly to the exact amplitude. Wallace points 
out the convergence of the series is worse for a Gaussian 
potential than for a Woods-Saxon potential. Gillespie et 
al. (GI75) have since examined numerically the Wallace cor- 


rections as applied to p - tHe scattering at 0.1 GeV, 
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0.5 GeV and 1.0 GeV using a central Gaussian optical poten- 
tial and find the convergence of the series to be excellent. 

In this thesis the concern is not with the relative 
merits of impact parameter formalisms, but rather in estab- 
lishing the accuracy of the Glauber eikonal approximation in 
reproducing observables from scattering by a nuclear optical 
potential. Several studies in the past have tested the valid- 
ity of the approximation as regards potential scattering, and 
so some justification is needed for this one. First some of 
the previous studies are reviewed here. 

Y. Hahn (HA70) has done Glauber eikonal and exact cal- 
culations on a non-singular Yukawa potential. Hahn consider- 
ed only one energy ( 1 GeV) and concluded that the approxima- 
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S.J. Wallace (WA70) whilst testing the now called "Wallace 
corrections" to the Glauber eikonal approximation, tested the 
unmodified approximation on Woods-Saxon, Gaussian and Yu- 
kawa real central potentials. 

F. W. Byron et al. (BY73) have examined the numerical 
accuracy of the Glauber eikonal approximation on potentials 
which are a linear combination of Yukawa potentials, and on 
a polarisation potential, and find Glauber's original cri- 
teria of validity unnecessarily restrictive. 

With the exception of Gillespie's, the aforementioned 


studies have been performed on potentials which have not 
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resembled nuclear optical potentials in the property of 
having absorptive or polarising effects. Gillespie's poten- 
tial had an imaginary term, but no spin-orbit term. A recent 
study by Brissaud et al. (BR75) rectifies these defects, by 
taking optical model parameters directly and also including 
the effects of the Coulomb interaction in the manner prescribed 
by Glauber (GL70), Brissaud et al find the interesting result 
that the approximation yields better results for the cross- 
section than for the polarisation. It will be seen later in 
this study that if one computes a polarisation transfer coef- 
ficient (P.1.C.) then this is reproduced less accurately 
Stik byetie: approxamatuona. Whe Pal. Cw used an’ sthis) study 

is the first Wolfenstein parameter, commonly given the nota- 
tion ae (RO67). Brissaud et al. alsofind the general rule 


of valididty 


a” « k/R 


holds qualitatively, but’ is unnecessarily restrictive for 
heavier nuclei. 

The main usage of Glauber theory is in the multiple 
scattering formalism, which is a combination of eikonal 
and phase additivity approximations. The eikonal phase from 
scattering off several centres of force turns out to be the 
sum of the single scattering phases from each centre, and so 
one can describe scattering from a nucleus using just the 
parameters of the nucleon-nucleon amplitude and the nuclear 


density. This is aesthetically pleasing since one can obtain 
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"macroscopic" results from microscopic input. To understand 
Whierestnie=eresults=irom the multiplesscattering are =valid, it 
is necessary to understand the validity of the impulse and 
eikonal approximations. The latter, it is hoped, will be 
clarified in this thesis. 

Tésts*havewbeen performed in the past ‘on the multiple 
scattering formalism; forexample on p - “He at. Various 
energies Auger et al. (AU76) find excellent agreement with 
the experimental data. 

Corrections to the multiple scattering formalism have 
been proposed and tested by Bleszynski et al. (BL76) who 
find a good fit, made excellent by their modifications, to 
the p- 58 Ns SCaULeG Ine Oatavate el eGeN.. 

The optical model can be tied in nicely with Glauber 
thoery by means of inversion of the Bessel and Abel transforms 
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tolpobtain an expression for V(r) (not unique) from the 
observed scattering amplitude. Brissaud et al. (BR75) doa 


calculation for 155 MeV p - 208 


Pb scattering where they use the 
multiple scattering formalism to compute the eikonal phases, 
and use these to obtain an optical potential. The potential 
So produced agrees well with the "best fit" Woods-Saxon 
potential, especially on the nuclear surface. A similar 


approach has been taken by Dymarz et al. (DY77), where multiple 


scattering formalism has been used to produce profile 
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functions for p - “he at 348 MeV, 650 MeV and 1050 MeV, and 
then instead of evaluating observables directly in the Glauber 
eikonal approximation, he evaluates an optical potential from 
which a phase shift calculation is performed to obtain the 
observables. The regular Glauber technique tends to produce 
rather sharp diffraction minima which are suppressed by 
Dymarz's method to give a better fit to the data. 


In this thesis we take an optical potential of the 


form 


2 
V(r) = (vy + iw,)(1 + pr’) oe 


2 
? 2 -ar 


The six free parameters are varied as a function of energy to 
Mimic the properties of the nuclear optical potential, in 
particular its strength, range and volume integral. This 
potential has the advantage of giving an analytic form for 
the amplitudes, but has the disadvantage of being unable to 
reproduce the scattering results from heavy nuclei; even for 
helium a Gaussian potential is not as good as a Woods-Saxon 
potential as shown by Frosh et al. (FR67), but nevertheless 
the results may be qualitatively fitted and so testing the 
approximation on these potentials is not entirely unrealistic. 

The potential contains a spin-orbit part, which there- 
fore compliments the earlier studies with central potentials, 
by having something to say about the spin-flip amplitude. 

Most tests before have been in the 1 GeV range, we 


concern ourselves with the “intermediate energy" range, from 
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0.1 GeV to 1.0 GeV in this study, and hope to shed some light 
on the successes and failures of the multiple scattering in 
this energy region. 

From spin-3 on spin-zero scattering one has four 
amplitudes. Accounting Rieis aa overall phase, one has three 
linearly independent quantities, this study appears to be 
the first to consider all three, taken here to be cross-section 
polarisation. and a polarisation transfer coefficient. 

The thesis starts with a recapitulation of the deriva- 
tion of the scattering amplitudes and observables in Chapter 
II and aslight modification is found for the spin-3 spin-flip 
amplitude. The modification does not always improve the 
accuracy and does not have a large effect on the amplitudes, 
formally, from the theory its presence is required and there- 
fore itis kept«-~inm Chapter Ili the results for a modified 
Gaussian potential are derived and also the associated Born 
results. Chapter IV summarises the numerical results and 
draws conclusions. Various results for cross-sections are 


given in the appendices. 
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CHAPTER II 
THE MATHEMATICAL BACKGROUND 
TO THE GLAUBER APPROXIMATION 


iba 


ii 


The Glauber approximation, like the W.K.B. approx- 
imation, is based on employing classical ideas into a 
quantum-mechanical framework. The Glauber approximation 
uses the classical concept of an impact parameter. The der- 
ivation of the approximation presented in this thesis unfor- 
tunately seems to mask this formalism. The essential point 
is that since we are using a semiclassical idea (i.e. we 
do not take just one impact parameter, we average over a 
probability distribution of them), we expect the approximation 


to be best when used under circumstances akin to classical. 


The "akin to classical" circumstances are that the 
flux is not changed much by the presence of the scatterer. 
This means that the potential must be weak in some sense, 
Tiiepattrcular: 

(a) The potential strength is weak compared to the 
incident energies. 

(b) The potential varies little over one incoming 


particle wavelength. 


This chapter. starts off by deriving briefly the Born 
approximation results for the scattering amplitudes of spin 
2 particles. The spin zero results are well known (R067). 
The reason for these Born calculations is that one may see in 
action the correspondence between Glauber and Born approx- 
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After these calculations the spin zero and spin 3 
ulauber approximation expressions are derived,: and a version 
for the spin-flip amplitude is obtained, which is slightly 


dyviterent trom that in the literature — —=“in particular (G59). 


Before performing any calculations, however, a recap- 
itulation on scattering theory may be in order, especially 
z 


for spin = scattering which is rarely touched on in simple 


terms in the standard quantum mechanics texts. 


Consider a particle of mass m and spin 3 being scattered 


Trom ampoceiural or the form 
V(r,0) 


where nos the position vector from the scattering centre to 
the incident particle and o& is the spin operator. Tite 


essential problem is, of course, solving the Schroedinger 


equation 
(974k “yatey = - zp ewlyite)), (1) 


The wavefunction for a spin 3 particle is represented by a 


Pu tie matrix, 


At large distances the wavefunction looks like a linear 
superposition of the incident plane wave and an outgoing 


spherical wave, 
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Wheresg\i; a'2 xe matrix called the scattering operator, is 

a function of © and g and is an operator acting on spin. 

a and b are two mambers satisfying a e4. b és 1 and describing 
the spin state of the particle. kK is the wavevector of the 


incident particle and k = |k.| since we assume the scattering 


Lemekacstice 


The differential equation (1) is best solved by turning 
it into an integral equation 


(Pats) = (3) eB + fot er eff pe o 


where G(r-r') is the free Green function which satisfies the 
equation 
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The solution to (4), which, when used in (3) satisfies the 
boundary condition (2) is given as 

eik|xr-r'| 

e(r-r') = —5 ———__.. (5) 


Thus, using this our equation (3) becomes 
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Assymptotically we have the expansion 
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Using this in (6) we have that at large distances from the 


Scalverlne centre, 
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where we have introduced the vector 


eee 
Comparing equations (2) and (8) one then obtains 
ai m= Sullen neg EN aS) 
M s)> . she “Sr Fvine Ye) ar. (9) 
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This is an exact expression for the scattering operator. 
We see from it that in order to calculate the scattering 
operator, we must know the wavefunction at all points where 


the potential is non-zero. 


Since any single particle scattering may be consid- 
ered to take place in a plane, it is worth at this stage 
adopting the convention to be used throughout the thesis, of 


fixing the y-axis to be perpendicular to the scattering plane. 


x 


Y, 


yh Ke 


y figures tL 
In this picture we have indicated, too, the choice of the 
z-axis, that is to be parallel to k,. Later we shall choose 


it to be midway between k, and ke. 
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THE BORN APPROXIMATION 


The Born approximation for the spin zero case is 
well known and well documented (R067), therefore the derivation 
is carried out here only for the spin $ case. The spin zero 
result may be obtained by considering only the diagonal ele- 
ments of the scattering operator to be the spin zero 
scattering amplitude. The actual wavefunction may, in theory, 
be obtained by iterating equation (6) until it converges. 
The Born approximation essentually uses just one iteration 


of (6) stanting from zeroy to eet 


Se) (0. : 
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We now must specify something about the potential. Let us 


assume it to be of the form 


V¥Cr,0) SMve(r) + Viet (rjo. Ly; (izs) 


which is used in optical model type calculations for nuclear 


scattering, and therefore of interest. m3] is then the sum 


ox two terms 


dr (13) 
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The first term is the scattering amplitude in the 
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where fig = n(k, - k,) = momentum transfer, and ee, is the 
zeroth order spherical Bessel function (AB64). The second 


term is evaluated as follows 
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THY GLAUBER-KTKONAL APPROXIMATION 


We study first the spin zero case, which is 
equivalent to ignoring the effects of spin, and then 
generalise afterwards. This derivation of the approximation 


follows that of Kamal (KA72). Starting from the Green funec-+ 
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The main difference between the first and second Green 
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have a bearing on the final result as we shall see later. 
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This is just the result derived by Glauber in 1959 
(GL59). The accuracy is improved by using the second Green 
function, which is where the propagator is expanded about 
the 2(ks + ke) direction. The derivation of the iscattenins 


amplitude is very much the same as above and is outlined 


here. 
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Note that this step is exact now, whereas before it was an 


approximation. 
For X(b) = X(b), we may do the angular integration, 
Since q.b = (k, = Kp) +b ‘ (72a) 
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This form is in agreement with the Born approxima- 


tion under eiX(b) __, 1+ ixX(bv) , and generally gives 
better agreement with the exact calculation. It is the form 
most used and is a manifestly time-reversal invariant result, 


since it treats kK and k symetrically. 
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We now go on to consider the case of particles with 


spin 3. 
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THE GLAUBER-EITKONAL APPROXIMATION 


FORS SP iINes sPROJECTILES 


The derivation presented here is the one which fol- 
lows from changing the axes in the manner that Glauber did 


originally, 4.46. 
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As before in the spin zero case this integral equation may 
be solved by differentiation and the boundary conditions decided 


by substituting the solutiom back into the integral equation 


yielding 
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We see that only the y term contributes to the second 
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Now, using the relation between the scattering operator M 


and the scattering amplitudes, we obtain 
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As Glauber points out in his 1959 lectures, the 
system of axes used earlier (k. parallel to z) is not the 
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(i) treats k, and k, differently, therefore making 
the result not manifestly time reversable from the outset, 

(ii) does not reduce to the Born approximation 
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In his article Glauber considers these points with 
the following prescription. 
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or (ii) since, if we assume en even parity potential along 
the z-axis, the integral vanishes. 

Hquatvor. (98) ls also nom loneerrexact bute thelconrec— 
tion term encountered by not making the approximation is of 
opeteiesa be i6l a relative to the remaining term, so we may neglect 
it to a good approximation. 
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we have to make the approximation q.r=q.b which requires 
k d 0° «1 . In the second we have to make bx Kk, = bk, 
which requires k d 0° «16; and rx k, = bx k, requiring 
k d 0 «1 only for even parity potentials. Since most potentials 
encountered are of even parity, then the second coordinate 
system appears to have the least restrictive approximations. 
(Both of course have the linearised Green function. ) 

Glauber points out the first improvement, which is that 
for equation (89) q.r=q+b exactly in this coordinate 


system, and then appears to say, since 
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This is adequate and it is true that the approximation so 
obtained is a better one than before, but inherent in 
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change coordinate systems at any stage one merely makes the 
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formalism, but in the above case there is no need for it. 
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factor cos 0/2. The scattering amplitudes resulting are 
co 


f((\e= ik] Jo(2kvsin 9/2) Pees b db and (120) 
BUC)s =—=1kcod 0/2 | 3, (2kesin @/2) feet) bed Gane (249) 


(2) 


One may argue that the presence of the cos 0/2 is splitting 
hairs, since we have already assumed small angle scattering 
and this may be true. Nevertheless there are reasons for 
keeping it, namely 

(i) it follows from the mathematics by starting with 
phewcecona- Green function and sticking to one coordinate 
system, 

(ii) it is a simple modification which does not lead 
LOeanyeturther jdiitveulty in) the calculations, 

(iii) it is necessary to make the spin-flip amplitude 
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reduce to that of the Born approximation as e —+1+ iX(b) 
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for small angles. 

(iv) it ensures that the polarisation vanishes at 180, 


even though one does not intend to press the validity of the 


approximation that far. 
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THE EXACT CALCULATION 


This is essentially aphaseshift calculation. The 
wavefunction is expanded in partial waves and the contributions 
from each partial wave toeach scattering amplitude are summed. 
This method of obtaining the scattering amplitudes is 
described rather well in (R047), a summary of which is given 
here. The exact results for the observables are required, 

Of course, to test the approximation. 


We have to solve the Schroedinger equation 
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The eigenfunctions will be in the "j" representation since 
the Hamiltonian contains a term coupling space and spin 
coordinates. The angular and spin parts cof the eigenfunctions 


are the generalised spherical harmonics, defined as 
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The boundary conditions require the wavefunction to be 
assymptotically a linear superposition of a plane wave and an 


outgoing spherical wave. 
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Using the plane wave expansion 
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Clearly demonstrated here is the difference to potential made 


by the spin-orbit potential for aligned and unaligned spins. 
These two equations may be integrated using the 


Numerovy (SH68) method and we may write 
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We may identify the scattering operator in here. Converting 


to the @-s basis, we have 
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The phase shifts may be found by matching the logarithmic 
derivatives of the numerically integrated radial wavefunc- 
tions to the logarithmic derivatives of the free particle 
wavefunctions at a sufficient distance from the scattering 


centre that the potentials have effectively vanished. 
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In this chapter the expressions for the scattering 
amplitudes derived in Chapter II are evaluated for a 
Gaussian potential of the form 
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venient expression for numerical caluculation of scattering 
amplitudes, a series in fact which converges reasonably 
quick vie 

The spin-orbit potential is chosen so as to have the 
Samemadecaye constant as the central potential.  Thissieere— 
strictive when it comes to fitting experimental data, but 
PULUSecdesinces1 tea SsHin dine with the philosophys of) thesanalytic 
form,of the’ spin-orbit potential being 1/r d/dr of the central 
potential -- it also makes the calculations much easier! 

The spin-orbit potential has also been taken as complex 
Since the optical model analyses have shown that this is the 
case for energies higher than 100 Mev. 

This potential in the Glauber approximation has been 
studied by Wallace (WA70), and his results have been reproduced. 
Wallace, however, considered neither the spin-orbit term nor 
complex potentials 
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The validity of the Glauber approximation as applied to just 
the central part of this potential has been reproduced numer- 
ically for the Glauber calculations. Included also in this 
chapter are the corresponding results for the scattering 
amplitudes evaluated in the Born approximation. The 
correspondence Glauber —-» Born as Boe ae SG ae BE7.€ sks 


noted. 


GLAUBER AMPLITUDE FOR A GAUSSIAN POTENTIAT 


We must first obtain expressions for the eikonal 


phases. The equations (93a) and (93b) give 
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Equations (102) and (103) now read 
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For the non-spin-flip part of the scattering operator 


equation (120) gives 
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Using uniform convergence we may integrate term by term to 
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The confluent hypergeometric function may be expressed as 
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The same procedure may be followed for the spin-flip 


part of the scattering operator. Equation (121) gives 
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If instead of a Gaussian potential a modified 


Gaussian potential had been used, 
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which are just the results obtained previously. 
If we use the Born approximation to calculate the 
scattering amplitudes for a modified Gaussian potential, we 


obtain the results 
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cos 0/2 missing from the 'g(0)' term of equation (196). 
The argument for the presence of the cos 6/2 factor 

multiplying the g(@) term gains weight here since then the 

first term in the Born series is reproduced exactly. 
Considering the unmodified Gaussian potential, the 


Born results are 
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Which again agree with the Born approximation except for the 


coss6/2 tactor in they \2(0)" term. 


OBSERVABLES 


Whereas from our calculations we obtain amplitudes f and 


g, experimentally we cannot measure these. Instead, using an 


initially unpolarised beam, all that we can measure from a 
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Single scattering experiment is the differential cross- 
section. With an initially polarised beam we may determine 
the polarisation from a single scattering experiment and 
the polarisation transfer coefficients may be determined 
from a double scattering experiment. 

If the Glauber approximation were to have the effect 
of altering the amplitudes by an overall phase factor, then 
the observables would not be affected. We therefore study the 
observables to determine the region of validity of the approx- 
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Differential Cross-Section 
This is defined as the number of particles scattered 


into a given unit angle per unit time for a unit incident 
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radial flux and using the Riemann-Lebesgue lemma 


lim joe ei K-£ dn = 0 (202) 
r= 00 
4a 


‘ie fas Sites - aes . 
ae res “he ns rennak seine 
; Ede) 


: endo a h SW naad & 


a ae 


bid Sask Preges githe aa aint, eu en 
eyo oo Fan athe Epi ak 


Sod ting bere “heat fy 


Parr mane -O7e0 = “a 0 emdksoraee % res ind T ot 
; VAN; - re a 2 va 
oy sity (iswve as ul mat igs’ a peeran ie : 


. rhea 

varie Sha euitewmag: aw Same: : 2 pl 2 home c+ idewrwode r BE P 
» = . - 7 

cs yao . = etm 

si he ys Bb cer: Yo -Mekgss 8aF »: iene ak cs colds ria 


» ang > oe 

ere) a SAZ saa a 

é i: —r 4 ‘ Pe a Met a. 4 a a & aia? 
janoos.p ADliLio Tes 2 Besant: She, 25 ror} ao -¥ oe 


i @ 


tor Ye 46) atla Ys TRE, rats = za 


Sal ae) metas Sit apideeed cahhenadeee OAR 
he he hy, 8 ant 


an | : j i 
: e 


. ery 
a 


2 \ —. : _ 
(00S) | lice + “7 


% 


ee 4ias) 


i] 


he Aelita Sr San 


ate i 


sen ib 


56 


(203) 
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The first term corresponds to the number of incident particles 
per unit area, the second to the number of scattered particles 
per unit solid angle. 

Thus, by definition, the differential cross-section 


is given by 
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This can be seen to be the "normalisation" for the scattering 


operator M in that 
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Polarisation 

The polarisation of a beam of particles is defined as 
the expectation value of the spin operator divided by the 
amplitude of the spin (3 in our case). Therefore 
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= 0 since all the Pauli matrices are traceless. 
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y is a unit vector in the y-direction, which is normal to 


the scattering plane. 


Polarvsation Iranster Coetiicient 

Phespolarvcavwon transter cociiicrents (Pa t.cas) 
give a measure of how the polarisation state of the incident 
beam contributes to the polarisation state of the outgoing 
beam. 
Mhemunderlvyines formalism of PolC.’s is describedsin (072) 
and (RO69). We will consider only the Ke coefficient, 
commonly called the first Wolfenstein parameter. The 
expression for this observable is 
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Therefore 
X= xcos O.— zZsin 6 (2172) 
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where the O*. refers to that component of the particle- 
Spin after scattering and o,. refers to the component be- 


fore scattering. 
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CHAPTER IV 
RESULTS 
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This chapter is mainly concerned with a comparison 
of the results obtained using the eikonal approximation with 
those obtained from an exact calculation. As mentioned 
in Chapter I, the energy region above 1 GeV has, for central 
potentials, been studied exhaustively, therefore this study 
is plimited =to the energy range from 0.1 GeV to 1.0) GeV: 

This so called "intermediate energy" range is witnessing a 
great deal of activity these days, and often the Glauber 
multiple scattering theory is used to describe the collisions. 

We choose here to study the elastic scattering of a 
neutral particle from a complex central and spin-orbit mod- 
ified Gaussian potential. Instead of picking arbitrary par- 
ameters for the study, parameters are chosen to reproduce as 
closely as possible the scattering of protons by actual 
nuclei, and then the charges are set to zero to simulate 
neutron-nucleus scattering. The nuclei chosen in this study 
are helium, calcium and lead, representing light, medium and 


heavy nuclei respectively. 


Choice of the potential parameters 


The potential parameters are chosen by taking the 
parameters obtained from matching a modified Gaussian potential 
to the "best fit" Woods-Saxon potential using the formulae of 
Appendix B, then using these as a starting point of a search 
for a potential to reproduce as closely as possible the elastic 
scattering data taken from the references in Table I. The 


resulting parameters are shown in Table II. 


; Laat oe ~*~. Ss vin . 

= a had » — - ¥) Lis earee Tae 

iat 6 ee a ee ea et 
a soy w aig TM nay pagal af 


a 
i : yh. Se - be : /* % a ae 7 — ry Sine 
ae 4 _ = = os F 
gota co ila eae ‘ieee @ 1: aul fh r Sey, ae 


yt a rat 


biasing: SA. cbt hetieLay tara. es *e at - ba 
levineo te? gel Vos ecid aciaes: yates sided Tost 
as. eisVen dnt totais ce “pat busta naa jaan 

(ed Gid.ad Ved Lo wey sgriex. Vaiaee She 8 bi ) 

) Beda ore dreooe. eynat “dais | teibedenrn e” vette 9 


wvselS of Walle Bas, ee -scecs ud Avie ad 14 inte ee 


= 


y wht ofttnesh 5¢, fie ae vane? bn inal elaine bel? “ 


Die 
fea teen Sigees ay chute oF oral: needa eW A 


a at 
ty ound hog Ppoorie Yeteier 2 ol aay ae nae 
' ta ne. e201 pie gost waive 2 
€ snl gute - Salioig To pas ant «Seer sal foe vai Bont 


> 


i. Lichen W ro te rut ohat tad) o A pal) 7 ay he 3 abe \d tee 


gqurow 3 ans 2eIceoe (OFF 254 0 oq mec: wo 


ores od item: edd sabe tala. gt yay foe 
; q ba 
aT fey orin of ‘3 a>’) omen ihe ah “sts sue oi setouart ee ne 


oe 


Sap. fel Lint « ae bt yaa riehae i er. Agods (uta Lae Lipo | ae bes 
Siererormme® to fons we 
‘ 2 - 


capa | hil eats, 9 uy . 
wh, erie ye neos a8 vit? eee Apitueay pale _ 
g/Paatou, caiaGwer HSrigtan « Ailirot ae mort benigcaa % sak 

rr on luswakl aie at Lauvulely ontag enteé-aboow ‘22. 
HoTeeh 2.20 2 eg.<gatarie » bu ed avid uae aad ¢ 
oiienia sar sldiiewy we cites ao wah oornies uw ia Samm 

| oft! oa edaal ia weenaesse palaps | mh 7 
Sina ie 


as 


iy 
a 


TABLE I 


The sources used for the experimental 
proton-nucleus elastic scattering data 


Target Energy(MeV) 


helium 100 
540 
580 
1050 


calcium 155 


155 
1044 


lead 1755s) 


155 
1044 


Observable 


cross-section 

polarisation 
cross-section 
cross-section 


cross-section 
polarisation 


cross-section 


cross-section 
polarisation 


cross-section 


Reference 


GO70 
GO68 
BO72 
AL75 
WI 68 


WI 68 
AL76 
WI 68 
WI 68 
BE 73 
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TABLE II 


The parameters used in the optical potential 


The potential used is 


2 
V(r,¢) = (votiw,)(ttper)e%F av euee tiw, | )(14pr*)e7%™ Ones 
Nucleus and V W V WwW 
re) fe) Sor SHor of [? 
Energy (MeV) (ev) (MeV) (MeV) (MeV) (fm) ( fm“ 
Helium 
100 VOnC By 15.0 5.0 iL BIOL QA 
380 10 80.0 @) 0) 0.0 ss (CEl 1.34 
580 473.0 sia e/S 0) 8.0 4.0 ik OM 1.34 
1000 -51.4 22785 Wen SS ol, lO hag 
Calcium 
155 9.0 14.0 et = OI On Oak 
500 =m 30.0 1.3 =(0) , 25 Orn Ont 
1044 -35.0 4O.0 ihe -0.4 Omi Ot 
Lead 
155 OeO fe) 546 ah ae) On OnOG Lin nOe 207 
500 O70 60 C Ono =O, 15 0.0644 0.207 
1044 7 O Z0=0 Orel 210) sak OOS 18) 207 
(1.0) (-0.5) 
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Fits were made for all the potentials used except 
for helium at 380 MeV, calcium and lead at 500 MeV. The 
parameters for these potentials were obtained by considera- 
tions of the behavior of the optical model parameters as 
a function of energy, and interpolating between the regions 
where experimental data had been fitted. Five fits are 
shown here (figs (4-7)). 

As expected from geometric considerations, the fits 
worsen as the ranges of the potentials increase, but in all 
cases the essential features of the curves, in particular 
the first diffraction minimum, are reproduced. Fits obtained 
with a more conventional Woods-Saxon potential would be more 
accurate (0OE76), but would lead to rather complicated expres- 
Sions for the scattering amplitudes, and so for the sake of 
simplicity are not used here. 

Generally the geometry of the potential was kept fixed 
for each nucleus, the exception being helium at 100 MeV where 
an increased range was required to reproduce the scattering 
data well. At 1044 MeV the lead potential obtained by 
minimising the chi-squared between theory and experiment gave 
a spin-orbit strength which was rather too low to produce a 
Mealvetic polarisation or P.1,C.,, the strength ol thevspin- 
orbit interaction was therefore increased slightly to the 


values shown in parentheses in Table II. 


Comparison of Glauber and exact calculations. 


For spin-3 on spin-zero scattering one has two 


complex amplitudes which, considering an overall phase acco 
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give rise to the three independent quantities. We therefore 
have to consider three observables in order to fully describe 
the system. We choose here the differential cross-section 


do : : : : 
AEN oO the polarisation P , and the first Wolfenstein parame ter 


xe HEU 
Ke (Zy Ne: 


The three observables are calculated here for three 
dpatercn tmucileil at various energies from.0.1 GeV tomo Ges 
and the graphs obtained giving exact and eikonal curves for 
them as a function of either the centre of mass scattering 
angle or the momentum transfer. 

What exactly constitutes a failure of the approx- 
imation is a subjective matter. One may be trying to repro- 
duce experimental data as in the multiple scattering use, in 
which case if the differences between exact and eikonal cal- 
culations are smaller than the error bars then one may say the 
approximation is valid. One may be wishing for a 5%, 10%, 

20% or "general features" reproduction, and so one should 

not conclude with a table of angles, energies and ranges 
showing decisively where one can and cannot use the approxima- 
GON: 

Despite the pessimistic attitude projected by the 
previous paragraph, a cursary glance at the curves will reveal 
that there is often a point where the approximation fails 
drastically (for example 80° for helium at 380 MeV (fig 8)). 
One may therefore draw the following Conclusions Ors uc 


various nuclei. 
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fig 4 


They tite tom thesexperimen ta iees pe “He elastic scat- 
tering cross-section data at 580 MeV and the fit to the 
experimental polarisation taken from 540 MeV, both eval- 


uated at 580 MeV. 


Helium 


580 Mev 
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fig 5 


The fit to the experimental ~p — ms 


ee elastic 


scattering cross-section at 1044 Mev. 


Calcium 


1044 Mev 
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The fit to the experimental p - 20855 elastic 


scattering cross-section data at 155 MeV. 
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The fit to the experimental p - Pb elastic 


scattering cross-section data at 1044 Mev. 
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Helium 

The cross-sections are reproduced rather well. At 
1050 MeV (fig 9) the approximation reproduces the first two 
minima, beginning to go astray at the third (see insert). At 
580 MeV the approximation may be said to fail at the first 
Manin recovering to fail sdrasticaliy sat 64°. For the lower 
energies of 100 MeV and 380 MeV the approximation does w ll 
(Uys We 80° on the cross-section. The polancisationvandwer a iac 


far leat or petore the first minimum for all enercies  nowever, 


Calcium 

The cross-section consistently fails just before the 
second minimum regardless of energy. The polarisation fails 
batdyeatetne; tirst dip at 150 MeV, sheer. Tt. C stall inves techie 
top of the first rise; corresponding to the same value of 
momentum transfer. For the other two energies the first dip 
Mie iew polarisation 1S Obtained, but the reproductionloiatne 


PeleCne Ls Noveappreciably improved. 


Lead 
Since experiments performed on lead, particularly 

at high energies are constrained to the foreward angles, it 

Uepoiminterest tomnove’ that in the case ol 1 Gey, stnes opin 

cation of the cross-section covers the entire angular range 

demanded of it. The approximation does rather well at 500 MeV, 


the polarisation and P.1T.C. being reproduced Significantly 


better than at 1 GeV. 
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As mentioned previously, to conclude with a table of 
angles, energies and ranges showing absolutely where one can 
use the approximation would be overambitious. However, if one 
decides upon a criterion of failure and then constructs a 
table, one may expound upon the relative performance of 
the approximation as regards these variables and hence no 
apologies need be offered for the presence of Table III from 
which one may draw the following general conclusions. 

With increase of energy, the angular range goes down, 
but the momentum transfer increases. Since the momentum trans- 
fer is physically more meaningful than the scattering angle, 


the approximation may be said to be betterat higher energies. 


For lighter nuclei the angular range is greater than 
for the other two, however the difference between the angular 


ranges of lead and calcium is not large. 


The cross-section is found to be reproduced to a 
greater angle than the polarisation. The reproductions of 
Polarisation anda P.T.G. are roughly the same, «tne P. 1. Cee being 


SlishtLy worse. 


Brissaud's (BR75) finding that the formula a° x k/R 
is unnecessarily restrictive for heavier nuclei is confirmed: 
lead fairs hardly worse than calcium. It is found here that 
the formula is also rather too restrictive for very light 


Nuclei eee nied Lum. 


Whilst increasing the energy from 500 MeV to 1044 Mev 
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TAB ieee 


The angles of failure of the approximation 


Nucleus Cross-section Polarisation Del aes, 
-1 -1 =41 
Energy ee q fm ON ens q fm ee q fm 
(MeV) 
Helium 
100 80 Pa S| 50 5 LO Lee 
380 80 4.5 42 PES 44d 245 
580 36 it 35 elf 55) Boe 
1050 80 Fi Lo ee 42 iis 
Calcium 
1555 32 iS 17 0.8 20 0.9 
500 20 ate, il) 1726 ae 
1044 2 hey, 10 age 5 8 sey 
Lead 
55 30 1.4 10 Om 8 0.4 
500 28 Phe dh Sail 125) aren 
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nduigeds Se) 
The differential cross-section for helium at 100 
MeV and at 380 MeV, calculated exactly and in the eikonal 
approximation. The solid line is the exact calculation, 


the dashed line is the eikonal approximation. 
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The differential cross-section for helium at 580 


MeV and at 1050 MeV, calculated exactly and in the ei- 


konal approximation. The solid line is the exact cal- 


culation, the dashed line is the eikonal approxima- 
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Diayshe 310) 
The polarisation for helium at 10C MeV and 380 
MeV, calculated exactly and in the eikonal approx- 
imation. The solid ine) is theyexactecalculadion, 
the dashed line is the eikonal approximation. The 


graphs are plotted as functions of momentum trans- 
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ses Ang 
The polarisation for helium at 580 MeV and 1050 
MeV, calculated exactly and in the eikonal approxima- 
tion. The solid line is the exact calculation, the 
dashed line is the eikonal approximation. The graphs 


are plotted as functions of momentum transfer. 
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gal Daeg 
The P.T.C. for helium at 100 MeV and 380 MeV, 
calculated exactly and in the eikonal approximation. 
The solid line is the exact calculation, the dashed 
line is the eikonal approximation. The graphs are 


plotted as functions of momentum transfer. 
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Presi 
The P.1T.C.. for helium at 5c0sMeVeandntO50esNev. 
calculated exactly and in the eikonal approximation. 
The solid line isthe exact¥calcula tion, thesdasned 
line is the eikonal approximation. The graphs are 


plotted as functions of momentum transfer. 
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the accuracy does not vary significantly; in several cases 

a greater momentum transfer is reproduced at 500 MeV than at 
1044 MeV (see for example figs 17, 18, 22, 23). This may be 
due to the fact that the potentials used try to mimic realistic 
optical potentials in which the real central potential changes 
sign somewhere below 500 MeV, also the reaction cross-section 
is smaller at 500 MeV, giving smaller imaginary potentials at 


Them _OWerenerey . 


The apparent early failure for helium at 580 MeV (fig 
9) occurs at the first diffraction minimum, however the ap- 
proximation is awry for only a few degrees and then recovers, 
failing again at a higher angle. At 380 MeV, where the 
potentials, especially the real central, are smaller, the 
first minimum is reproduced more accurately (fig 8) and so 
the variance between exact and eikonal curves is not consid- 


ered (by the author) a failure. 
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fig 14 
The differential cross-section for calcium at 
155 MeV, calculated exactly and in the eikonal approx- 
imation. The solid line is the exact calculation, the 


dashed line is the eikonal approximation. 


exact 


reps ft -eikonal 


Calcium 


155 Mev 


oH 


: ie Saal 4 aed ¥ “as 
dy? Upege ang ue hea —— 


: ; 
PP 
ai) * : ' , wos 


a 
areas bar tee 


An 


<2? % 
a 
i a 


gO 


aT A 
«* = 
mS 
s & _ au “ - ra i? wd ~ -< . , 
ts MRINisg Hk ng ifese-euato igignesers if add 
} 


- ‘ , ———— 2 an ; >’ 
Pe: i ee 7 Dal De one a Oe ViFo Bee » ae f a Sig m4 vex 00 a 
os ke: eer eh a sat th ae iS oe ae 


pala iss eat wip al agli Oio4 oaP opie eee | 


rewt 


_— = pi.Lenal 


neaWel- RS) 
The differential cross-section for calcium at 
500 Mev, calculated exactly and in the eikonal 
approximation. The solid line is the exact calculation 


the dashed line is the eikonal approximation. 


exact 


eikonal 


Calcium 


500 Mev 


ao 


SA 


ie 
~~ 
- 

~~ 


7 | 4 \ ! ; 
1 7 \ . i <J)" 
4ats- [so tot notfose-nastry Lottossettib eat 


| v1 1 

_ ' ‘ , i 

“KOT s Taro le ovr of Ove yitorxs Sevaltofss yrot do) | 

& ? “ 7 . 
$901) fuote» ¢ Jax? os of 


ee 


“4 jrottamtxotgqs lenotils eit ef stil. bedash 


ote em ss 5 £ I Og i Pe z «28 i? yn i 


= ¢ — 
| 


eaa¢? 


atas -=7= + sel None) 


cha ako) 
the diiteréntial cross-sections or calovumes | 
1044 Mev, calculated exactly and in the eikonal approx- 
imation. The solid dine is-the exact. calculation. ge 


dashed line is the eikonal approximation. 


22 


= 
= 
ra 
3) 
qo 
w 
oO 


1044 Mev 


aa ee CK Ona 


¢ - ar aD i 


te! 


~ * i - 
oad AR ah ie a | 


12221 : : 7 
: : jr 


[nnoRs a= += = = 


: 50° Sey 
: 
apem G2i te wloslao tot nofs KeLISLOg BAT 
\Cronz6 Ostvatwolas vol MOL brea 
Bivioiso toaxn) ery el eanis biiow onl ..noliemixougs 


ee 


atl HO LF at KON t lecodie off 2! ull Botesb eff 


ardtenert act coon to Solfoeavwl 3 es be7Tsoig STs Stiga 


0.4 - 
"7 


ne aWed Gk y/ 

The polarisation for calcuim at 155 Mev, 500 Mev 
and 1044 Mev, calculated exactly and in the eikonal 
approximation. The solid line is the exact calculation, 
the dashed line is the eikonal approximation. The 


graphs are plotted as a function of momentum transfer. 


150 Mev 


500 MeV 


exact 


eikonal 


1044 Mev 


Te 


7 eae a ad . - 
Ee 


¥ 


_— 


:. = eas oe 


7 


ad 
<A. i 
‘ 


, . z von owe vou cc f ye it: mt ft 
» 7 ‘ 


& 


: ; ro -~_ + *. ” 
Fay at ent ia elena aia ‘5 
7 i “ 


- ® 
= , telornet mut anos to 
| : 
, d 1 | 
tf 


\) 
— 


mie ,nolveivsifo foeth. tit 2! 
S| § 


“arodte ed Bs | 


ay 


a7 | Pe] tO 2 vou 


} i’ 
otgi Biloa 


/ ‘> : 


ae 
bol ee 5 


i 


f ont 


’ rr. 
‘eOurge Lenore s/t oi oom yltoaxs peteivoles Wem eos 
! ¥ . 


? 


; 

} 

f-\ | 
aT  teoiven 
4 


iS end C\ ophas 4 


m ifonut\e/as Bebe btr vas ! 


Lp a 


155 Sev 


iwiee ake: 
The PeT.Ce) for calciumeat 155) Nev, 5S00RNeveand 
1044 Mev calculated exactly and in the eikonal approx- 
imation. The solid’ line ve the exact calculation jeuae 
dashed line is the eikonal approximation. The graphs 


are plotted as a function of momentum transfer. 


99 


1044 Mev 


500 Mev 


exact 


eikonal 


i 
a8 i. a Der 
A 


Lininoel.£ ws — ee - & 2 S= 


Gi alt 
, | wi ahteett meoto Laiinetertip odT 
: rat Q . Lanois afy mt bas yifeaexs betetvoisn vom 
Defeah 0 uno rau ao foune ot af omtl bifos edt 
i * ren benoshe edt al 


eo ia 
} 


a 
—_ 


— 


a 
ee ae ee 


——— 
= 
~ 
4 7 


sine IS) 
The differential cross-section for lead at 155 
Mev, calculated exactly and in the eikonal approximation. 
The solid line is the exact calculation, the dashed line 


is the eikonal approximation. 


Od 


Lead 


155 Mev 


lyr@ 


' oer 8 hf . OS sit 
: a ' 
* 20, +s bbe 10% nbstoGm- 22019 Lsitasierilh att 


7 7 uJ . 
30 sd Kenoirr edt oi Bae “(fosxe betaluolbo .vek 50 


mB 

‘e suit biloe oft -.qofrandt 

- ae Jnoit andronege | lxwAre eff ef anil betieeb 
{ \ 


| 


fe ,~NOltelvoiss Tacxe2 st’ 
: eum | 


a 
‘J 
| 


i] 
i 
i . | 


Vv 


Wess (40) 
The differential cross-section for lead at 
500 Mev, calculated exactly and in the eikonal approx- 
imation. The solid line is the exact ‘calculavion, the 


dashed line is the eikonal approximation. 


103 


Lead 


500 Mev 


exact 


eikonal 


inp 07 


ape Wy. ed Greer) 


jeteauiLews dour i 


} 
fs ait! | 


fi wot noltoar-asers trevaat 


y/ A 
Ga ferouie adit ai bers Qitvete, bg 
MG ssRLoiao Jowxe wit ef oalli 


ee 
iemoivamizowpye Lancdle 


ne rhied  rell 
The differential cross-section for lead at 1044 
Mev, calculated exactly and in the eikonal approxima-. | 
tion. The solid line is the exact calculation, the 


dashed line is the eikonal approximation. 


105 


1044 Mev 


at 
fay) 
S 
fe) 
wd 
o 

! 

l 

1 

I 


inal 


| 
Vat EQ? 


, ' . 4 | 
’ ' 
F .o ‘ 
ogy Yo a 
‘el eae 
} if 
f 
fl Ht 
; 
Bi ' + 
i 
P04 oi fl 
i 
— a | ab 
‘ a , 4 f } 
} vr és | if { b-¢ ii 
7 A en Meee | * 4 
= ' i ’ 1 . 
. us 
A ¥e% ie a as rah, be besttf 191 cab gt ae toda ty act \ 


U ‘f ; 


paaTs we taht ‘aad nz Bes vfeeuange here Sn iSeo vis alg 

“neds htve tne tence sit al oes biibg itt Samat 

-fegaess ed? |! .mot} solno1qa2 Lononts aie! a ls emit hiviews 
| ‘ 


‘lanes metaewoa to woos é bet tolg sis 


=~ a= » » 


ae 3 oe 


t ' 


D 2% 


: nay + 
ow aT 


Tees 
The polarisation for lead at 155 Mev, 500 Mev and 
1044 Mev, calculated exactly and in the eikonal approx- 
imation. Thesolid line is the exact calculation, the 
dashed line is the eikonal approximation. fhe graphs 


are plotted as a function of momentum transfer. 


an0 7 


> 

o 
= 
= 
= 
fo) 
a 


500 MeV 


eikonal 


m~ Se — — HS |S 


, | v 
ss Av 
Yoo ove a : . 
, 4 5 
- y : 
ee ee . mona orn 
° | 


© : ™~s 
2 vad OF a2 
; heteabte. 4 ‘ 


ene ee Waele § kis ult widlands 5 
aan ey ..40/ arkxotggs | Jsandie 


PTT ideale sre Lo sade £ es 


7 Ve 


cBeh 725 
The P.7.C. for lead at 155 Mev, 500 Mev and 1044 
Mev, calculated exactly and in the eikonal approximation. 
The solid line is the exact calculation, the dashed line 
is the eikonal approximation. The’ graphs are plotted 


as a function of momentum transfer. 


109 


500 Mev 


— = 
seeded) i 
—-=—-. 


o-- 
- 
-- 


-0.5 


eikonal 


— eo 


> 

© 
= 
=F 
=F 
co) 
ban 


ce ct Mg ili a Re AE —— 
: : 


a oS ae 


; 7, f | Mids: = 
<4 ¥ tis ie, uy qos pte Me 


phn 4 gaa ket ; shila _ 


bn iin! nae | 


110 


Comment on the amplitudes 

It is observed (see figs 24 - 27) that in the 
replication of an amplitude by the Glauber eikonal approxima- 
tion, the relative error is usually greater for the smaller 
of the real and imaginary parts. This may be explained in 
terms of the Born-series analogy; if the actual error in each 
term is about the same, then the relative error in the result 
will of course be larger if the series converges, through 
cancellations, to a relatively small value. 

The effect of the greater relative error in the small 
part of the amplitude explains why the cross-section is repro- 
ducedsmore saccurately than the polarisation or) the Pat.ce, 


which are sensitive to the individual phase of each amplitude. 
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Tne et#ect of this factor is seen to be quite small. 
In the example given (fig 29) the factor is seen to improve 
the reproduction of the polarisation, but worsen the reproduc- 
TvOnmOtetiesP wr. c.. <n Other cases stestedea Similars scomilor 
"random" amelioration/deterioration was observed, making it 
difficult to say anything about a practical reason for retain- 


imewor dvecarding the factor. 
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The exact and eikonal non-spin-flip amplitudes 
tabulated for helium at 100 Mev as.a function of both 
centre of mass scattering angle and momentum transfer. 
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tabulated for helium at 100 Mev as a function of both 
centre of mass scattering angle and momentum transfer. 
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Table VI 
The exact and eikonal non-spin-flip amplitudes 
tabulated for Noite at 500 Mev as a function of both 
centre of mass scattering angle and momentum transfer. 


The units are fm. 


taken from calcium at 500 Mev 
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(The exponents of the numbers are shown in parentheses.) 
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Table VIl 
The exact and eikonal spin-flip amplitudes 
tabulated forscalcium at 500) Mevias astunction of, bow 


centre of mass scattering angle and momentum transfer. 
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SPIN-FLIP AMPLITUDES 


taken from calcium at 500 Mev 


TABLE OF EXACT AND ETKONAL 


(The exponents of the numbers are shown in parentheses.) 


Few WwWNhH DH FP PO OO OO 


Re(g) 
exact 


ORO 


.027(0) 

.548(0) 

-575(-2) 
.667(-2) 
.057(-2) 
.186(-3) 
.608(-3) 
»284(-3) 
-335(-3) 
-7. 


672(-4) 


Re(g) 
Glauber 


Onc 
4, 344(0) 
1.465(0) 

-1.673(-1) 

-6.769(-2) 
1.759(-2) 
330g) 

=i Jabal 5h) 

-1.175(-4) 
8.691(-5) 
2.871(-6) 


Im(g) 


exact 


0.0 
9. 
5). 
8.724(-1) 


174(0) 
035(0) 


1.054(-1) 
Delole-2) 


.220(-2) 


1.263(-2) 


.117(-2) 
-890(-3) 
SEV ACS") 


Im(g) 


Glauber 


00 
A47(0) 
2p To 
P30 (ei) 
»259(-1) 
Cones) 
Po (ea) 
.785(-4) 
Poct (es) 
HW (6) 
.872(-5) 
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Dige2s 
The polarisationsand| P.T.C. curves, calculated 
exactly and in the eikonal approximation, showing the 
effect of the cos 9/2 factor. The calculations are 


made for helium at 580 Mev. 
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CHAPTER V 


CONCLUSION 


aal 


The Glauber eikonal approximation has been studied at 
intermediate energies i.e. (0.1 Gev to 1.0 Gev). The study 
has been undertaken using a modified Gaussian Otel as: 
rather than the more conventional Woods-Saxon potential, The 
Porumerwis chosen for sits simplicity. 

Most earlier analyses were done using only central 
real potentials. In this study we use a more realistic 
complex potential with the added important feature of includ- 
ing the reather inevitable spin-orbit interaction. This 
enables us to stucy the behaviour of spin dependai.i obser- 
vables such as polarisation and polarisation transfer 
coefficients. From the results of the study we may tentativ- 
ely draw some conclusions. 

It seems that the performance of the approximation 
does not significantly improve as one increases the energy 
from the 500 Mev range to the 1000 Mev range. The optical 
potentials used in this study mimic realistic optical poten- 
tials in being weaker at the lower energy, this may explain 
the Low energy success. 

Generally one finds that the approximation fails at 
2 smaller momentum transfer in the 100 Mev region than for 
the higher energy regions studied. 

On average the polarisation and polarisation transfer 
coefficient reproduction were much the same, the approxima- 
tion tending to break down at momentum transfer very roughly 


around 60% of that at which the reproduction of the cross- 
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section fails, but this varies congiderably from case to 
case. 

It bears repeating at this point that what constitutes 
a failing of the approximation is a subjective matter. 
Nevertheless a table where the significant deviations set in 
is given in Chapter IV (table III). From the table he 
Following trends emerge. “itis seenrthat for a specitic 
nucleus, increasing the energy tends to decrease the angular 
range of the approximation, but increasesthe amount of 
momentum transfer reproduced. The lowest points of failure 
Hounduraethicestudvyeareate10 wand=cs 5 fm + for the centre 
of mass scattering angle and momentum transfer respectively. 

An extra factor was found multiplying the spin-flip 
amplitude, which was not present in previous works. Formally 
the presence of this factor (cos 0/2) is desirable since it 
leads to the expected reduction from Glauber to Born 
approximations in the limit of very weak potential strengths 
andwa Vanishing polarvea tuonwal eloUs ss ines eltecG, ote une 
Paccor= One tne accuracy of the approximation was found to be 
small, it's presence neither consistently worsening nor 
improving the approximation. 

For scattering off helium the results of (DY77) that 
the eikonal approximation produces diffraction minima which 
are too deep are reproduced. Their study also uses ~The 
profile functions from the multiple scattering to obtain an 
optical potential on which an exact calculation gives a good 


fit to the data. This indicates that the simple eikonal 
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approximation is not doing too well in the intermediate 
energy region for helium. However as shown in (GI74), the 
eikonal approximation does well in this region with the 
Wallace icorrecitions. 

For calcium this study has shown one may trust the 
eikonal approximation only as far as the first minimum if 
one considers polarisation, but the cross-section alone is 
generally reproduced up to the second minimum. 

For lead the approximation seems to be doing 
surprisingly well, reproducing four minima of cross-section, 
polarisation and polarisation transfer coefficient fail early 
for the lower energy but are reproduced well for the higher 
energies -- spectacularly so at 500 Mev. 

Thus we may conclude that the Glauber eikonal approx- 
imation must be used with caution at intermediate energies, 
especially where the observable in question fluctuates rapidly 
Since the approximation has difficulty in reproducing minima. 

Whether one should use the simple eikonal approxima- 
tion or use corrections to it for a certain case is an open 
question. The accuracy in the case of helium is significant- 
ly improved by the Wallace corrections or by means of an 
intermediate optical potential, however much simplicitly is 
lost; and so the final decision must be made by the require- 


ments of accuracy imposed upon the approximation. 
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APPENDIX A 


Total Cross Sections 

ise am ataits appendix are calculated the total cross 
section, the scattering cross section and the reaction 
cross section. The physical interpretation of these quartities 
is well known and is expounded upon in great length in almost 
any book on quantum theory -- for example (R067). 

A brief summary of the quantities is given. A particle 
is said to have been scattered during a process, if its state 
after the process is different from that before the process. 
Here the word state refers to both the internal state of the 
particle and also the momentum state of the centre of mass of 
the particle. A particle is said to have been elasitically 
scattered ifit has been scattered and its energy is the same 
in the final state as in the initial state. 

The total cross section is defined as the ratio of the 
number of particles scattered per target particle to the 
number of incident particles perunit area. Itisdenotedby , 

The elastic scattering cross section is defined as the 
number of particles elastically scattered per target DaLLicLe 
to the number of incident particles per unit area. It is 
denoted by Osc 

The reaction cross section is defined as the ratio of 
the number of particles scattered, by processes other than 
elastic scattering, per nucleus, to the number of incident 


particles per unit area. It is denoted by Op 
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The above definitions imply the following relationship 


On =OR +05 (226) 


The total cross section gives a measure of the number 
of particles deflected from their foreward motion and is 
therefore related to the foreward scattering amplitude. The 
relationship is the Optical Theorem (R067). 

4tr Im 
On SS eee ie 100) (227) 
k(2s+1) 
where k is the wavevector of the particle and sis the 
spinor thedparticle. 
Ltr 
Ong=w=— sim tO). (228) 
eee 

The scattering cross section is obtained by summing 
the probability of scattering into a certain angle over all 
angles. 

The probability density of elastic scattering at 


angle @ from spin state » to Spin state v is given by 


| m ,(@)| ° (229) 


Thus to get the total probability of elastic scattering, we 


must average over the initial states and sum over the final 


states. 


. 
| Is t001? da (230) 
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Q 
| 


sc = | { [sqey hens (ig(o)|* + j-ig(e)|* + Ir(0)|?} aa (231) 
f' |z(0)|* + Je(@)]*) an (232) 


H} 


Where 200) 1Samsin Ode dg . 


Oo. is most easily calculated by subtracting the above ex- 


pression for Os from On : 


The exact calculation of cross sections 
We obtain the total cross section from unitarity. 


Specifically using the Optical Theorem 


Op = a Im £(0) (233) 


on our amplitudes 


(234) 
(tO. 1 Y aries oes jy + 1619) sin Ma ee cos Q) 
‘ey 
°° es 
g(0) = erify - 07153 yo! (cos @) (235) 
1=0 
which gives 
Ke 3) 
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It must be remembered that in optical model analyses 
a complex potential is employed, which gives rise to phase 
shifts which in turn will be complex. 


To obtain the scattering cross-section we use equa- 


(aes) 
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expression for f(0) and g(@), we get 
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e2idy - 92s = 2ife sind, ~ eiby sin 1 } (243) 
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The exvression for % may be simplified, getting 
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The Glauber Approximation 


The scattering cross-section 


We evaluate the expression 
Z. Zz 
if [f(o) |“ + |e(e)]|~) an 


in the eikonal approximation. 


Evaluating the expression 
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The term in the brackets may be evaluated using an approxima- 
tion. The integration is of course carried out over the angle 
made by the final wave-vector, which means kK, is effectively 


a constant for this integration. The brackets become 


f (252) 


We are integrating over the sphere §S at present, but for 
small angle scattering integration over S is approximately 
the same as integration over the plane P. Since in the 
eikonal approximation we always assume small angles, then we 


may make the replacement, obtaining, for the bracketed term; 


oi (BBY) [oriite: (BB) —. (253) 


(The volume element in the plane is a“K., on the sphere it 


is |k,|“an,.) 
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We can simplify to obtain 


_+(b-b') 2 
(2m)* S(b-p') = “LD (-p"). (254) 
kK 


Substituting this back into (251) gives 
fea Zero 
ne [ir doa (255) 
Evaluating the expression 
of = { leconl? ane (256) 


using the expression for g 


g(@) = - # [eta sing [,(b) 4°, (257) 


which comes from equation (111), we get 
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Now, using the same approximationas before for the 


term in the brackets, we get 
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The optical theorem (228) gives 
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Adding (255) and (261) we obtain 


Ose = [Real + | Pi(b)|? sintg ) a2 . (2644) 


The reaction cross-section is obtained by subtracting equa- 


tion (264) from (263) to obtain 


O, =| (2ReM(b) - | CD)? - |h(w)]? sin? ) a2 (265) 
={ ( (lb) + Px(b) - Py(bPe(b) - | ,(b)]* sin®g ) oe 
= [(1-]1 - Myce]? - [rycel? sin®g | a’ . (267) 


Equations (263), (264), and (265) give the eikonal results 
for cross-sections. They may each be evaluated in a similar 
manner, and so only the result for Op is obtained here. 


The other two results are quoted. 


Reaction Cross-Section 
Making the assumption of azimuthal symmetry for the 
potential in equation (267), and using equations (102) and 


(103), we obtain 
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Now, using equations (155a) and (155b) this becomes 


Z 
eae Laneige 
ce {Jeo (bA .e av ) | 


| 
ll 
nN 
x 
co) 53, 
TES 
I 
oO 
& 
QO 
4) 
I 
& 


(270) 
(a Wg 
+ % |sin (paje"* )| } J: db 


vad a iz 

-ab 49) ’2 2 

an |v ab [1 — eA o& eAGe {cos (ba e7 ™ )eos (ba%e =P ) 
: Ve: 


2 2 
+ dsin (bale * )sin (paxe™ % ) 


9 
I 


- ob* - «b* 
% 
> - elA tax )e {gs (b(A +A )e”™ ) 


R nf db 


AOL OE Aare Lcpea Ma “be 
+ cos (b(A-A* )e ) +$cos (b(A,-AX )e ) - cos (b(A +A Je 


(272) 


Introduce 
of =) AIG) Gal heel Cr Vale (279) 
ip [(A_+A*)+ib(A_-A*)] “wD 
ORC Aa) = 5 |» (2 = Ze Cn Cuben a’. Sues 
(274) 
_ ALA tag )+ib(A tag )) 0% 
i = Rated. ax ))T ~urb* 
= 1b db S [omer l (act ag rin(a,- ee 
“ a G275)) 


2 
; r_-arb 
~ pr [(A tat +ib(A+Aa%)] “e 
Since the integral couldn't have been done analytically, we 
had to expand. The series is uniformly convergent OTe t 


0 < bgeo, thus we may interchange sum and integral to obtain 


te aout tat Bi 
Seg | st. 
(s¥s) | 7 leone oe 


ats. 
ages asa lta wna coat ( “ge a ee 


oe DPT AYRE T. . ass 
(eye) hhh Ali 
om argorgmang shee its. 
cans) peek a 
se ora ah AR o e sea 


fe ee on aaa a 


137 


an expression which still defies analytic integration, how- 
ever upon further expansion the integrals may be performed 


to obtain the expression for Op: 


[E] 3 (276) 
eo. pli (AetA* ) 7 5 °5(-7) Sa) 
‘ vee Sa earrer ss Ss | 


The total and scattering cross-sections may be evaluated in 


a Similar manner to obtain the results 
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APPENDIX B 


Matching Woods-Saxon Potentials to Gaussian Potentials 

We attempt here to give a prescription whereby, given 
a Woods-Saxon potential,we match it with a Gaussian potential 
so that it will have similar scattering properties. 

There are several physical characteristics that one 
may wish to equate. They are 
(i) the volume integral, 
(ise uLhe fr. nas. mradime, 
(iii) the potential gradiants at the nuclear surface 
(iv) the depth, 
(v) the strength at the nuclear surface. 


The analytic forms of the potentials matched are: 


(v_+ iw ae + pr*) + (Vv + iw yen ®t (4 + pr’) (ens Jt 
O O = Dee SuOn S.O% =— = 
and 
(v_+ iw_) al 
O O ids 
1+ broke a an, i iy ara) ar ee EEE ik - 


Often only expressions for P= O will be given;however, all 


formulae that are used in this study are given. 
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Modified Gaussian 


ee sole Tink -1 
(r5 Airy ia Fetex ey) ana 


Cine) Potential Gradiants at the Nuclear Surface 
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= Tq(pR(4a - R) - 1) 


r=R 
Gaussian 
2 Zz 
d -a“r bes -«R 
ae ve | = -2«Rve 
r=R 
Modified Gaussian 
2 2 
d 2 -ar _ -R 2 
APM cme) 2 | = -2ve RCQK=+2 (Re 818) 
r=R 


(iv) and (v) The Depths and Strengths at Nuclear Surfaces 


These are found by evaluating the potential function 


at r= 0 and r= R respectively. 


General 

It was found during the study that matching potentials 
Was Satistactary. for light nuclei. The criteria for a Woods- 
Saxon potential to be well matched by a (modified) Gaussian 
is that R/a should be small. The critical value of this 


parameter is around 6 above which the match is unsatisfac- 
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tory, For example, consider the best matches obtained for 


Netume(R/a = 3.008 and lead CR/ ates 
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Ti eue9 
The sucesses and failures of matching short and 
long range Woods-Saxon potentials by a modified Gaussian 
potential. The solid and dashed lines are the modified 


Gaussian and Woods-Saxon potentials respectively. 
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APPENDIX C 
Summary of Tests Performed on the Calculation and Program 
Three separate checks were performed, each to calculate 
the scattering amplitudes from a different potential. The 


potentials were 


% 2 
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STO S70. 
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(atalau)) Vvoe (1 + pr”) 


These give rise to expressions for five scattering amplitudes: 
£,(0), £2(0), £,(0), €,(0), and g5(@). 


The following reductions are found: 
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Three separate Born approximations were performed, and 
the results compared to the first term (r= 1) in the series 
Oca ri 


ae &4) and Bo : Providing one keeps’ the 


ie Byte 36 
cos 0/2 factor in the g terms, then it was found that the 
r= i term was identical to the first term in the Born series. 
In his thesis (WA70) Wallace gives the expression for 


the scattering amplitude from a potential of the form 
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This expression is reproduced by; 
ip iN Ag 8 1) Aer st e 
By oS One mesons Bric Ccde feet 


fr tp 7-0 


From each of the three separate calculations 
three computer programs were constructed, and tests were 
performed upom the numerical results obtained from them. 

Under the reductions described in the previous 
section of this appendix, the programs produced 
identical results. 

For each of the programs the first (and simplest) 
term was evaluated. This number tallied exactly with 
the results obtained from a hand calculation using the 
first Born approximation. The calculation was repeated 
for angles from 0° to 180 in steps of 5° and the tally 
between results was still found to be exact. 

For the first program the third and fifth terms 
were evaluated by hand, these, too, tallied with those 
obtained from the program. 

If one chooses to match a Gaussian potential by 
a Woods-Saxon potential, this may be done very closely 
(as opposed to the other way round!) and upon using the 
sectton of the first program to compute exact scattering 


amplitudes, then the amplitudes were found to be close. 
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It was found that in the high energy limit, and 
for small angles, that the Born and eikonal results were 
close, as they should be. 

Wallace's (WA72) results of comparison of exact and 
eikonal calculations appear to be reproduced, although 
Wallace only gives a graph, there is no visible deviation 
between the two sets of results. 

Y. Hahn in a series of three papers (HA69) (HA70) 
(HA73) has done fairly extensive tests of the eikonal 
approximation, and has numerical results tabulated for the 


scattering amplitude obtained from a potential of the form 
2 
-“r Z 
Vac ie Cee) 


His eikonal approximation results are reproduced by 
programs (ii) and (iii), his exact calculations are almost 
reproduced, the small difference between exact calculations 
being thought to stem from a small error in his calculation, 
possible (but not checked) taking too few phase shifts. 

So far the tests have always been on the non-spin-flip 
scattering amplitude. Noreferences found seem to give numer- 
ical results for the spin-flip amplitudes, except (ER75) 
which uses a Woods-Saxon potential. 


We discuss next a test performed to check the spin- 


Clips amplitude. 
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Relation between the two scattering amplitudes 


For a spherically symmetric potential 
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A\(o) is the eikonal phase function, defined by 
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Therefore 


ae fe 
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Cc s 


Note that if one has used the definition of the e2(0) with 
the cos 0/2 factor, these formulae would take a more natural, 


Simpler form 
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This is yet another argument (albeit weak) for the presence 
of the cos 9/2 factor. 

These relations were tested numerically and found to 
holdee. Thissis comforting, butssnows only thatec(O)eis being 
calculated correctly within an additive function of all the 


variables except ve and we! 
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